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Abstract 


Hypersonic flow over spherical dome protuber- 
ances was investigated to determine increased 
pressure and heating loads to the surface. The 
configuration was mathematically modeled in a time- 
dependant three-dimensional analysis of the 
conservation of mass, momentum (Navier-Stokes), and 
energy equations. A boundary mapping technique was 
used to obtain a rectangular parallelepiped compu- 
tational domain, and a MacCormack explicit time- 
split predictor-corrector finite difference algo- 
rithm was used to obtain solutions. Results show 
local pressures and heating rates for domes one- 
half, one, and two boundary layer thicknesses high 
were increased by factors on the order of 1.4, 2, 
and 6, respectively. However, because lee-side 
pressure and thermal loads were reduced the two 
lower height domes did not experience any net 
increase in total loads. The total loads on the 
higher dome were Increased by twenty-five percent. 
Flow over the lower dome was everywhere attached 
while flow over the intermediate dome had small 
windward and leeside separations. The higher dome 
had an unsteady windward separation region and a 
large leeside separation region. Trailing vortices 
form on all domes with intensity increasing wit;) 
dome height. Discussions of applying the results 
to a thermally bowed thermal protection system are 
presented. 
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viscosity (eq. 6) 
density 

stress tensor (eq. 2) 
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Introduction 


Protuberances in the form of large-diameter 
small-height spherical domes often occur on hyper- 
sonic vehicles as a result of hardware design or 
thermal bowing* Hypersonic flow over these protu- 
berances can produce increases In local wall 
pressures and heating rates. Interest in this flow 
configuration and computation of the pressure and 
heating rate distributions on them has increased 
because of proposed alternative metallic Thermal 
Protection Systems (TPS) to replace the fragile 
Reusable Surface Insulation (RSI) on the original 
Space Shuttle. 1 * 2 * 3 Titanium and/or superalloy 
designs that are weight competitive with RSI and 
offer the advantages of mechanical fastening, 
longer life, and increased damage resistance have 
been fabricated. 


e 

h 

K 

k 


L? .t-rj ,k 

P 


r 


o 


T 

t 

u, 

7 


v, 


w 


X* y, 2 

5 fp 


e 

n, K 


NomencI ature 


specific heat at constant volume 
(eq. 5) 

internal energy 

dome height above flat plate 

thermal conductivity of air 

variable grid concentration 

parameter (eq. 14) 

differential operators (eq. 19) 

pressure 

gas constant (eq. 4) 
radius measured from dome 
centerl ine 

radius of flat-plate/dome 
i ntersection 
temperature 
time 

velocity components 
functional relationship between 
the physical and computational 
domains (eqs. 13 and 15) 
coordinate directions 
flat plate boundary-layer 
thickness at x - x f 
one-half the fillet width 
transformed coordinates 


However, the use of metallic TPS can introduce 
some physical phenomena not present in the RSI 
system. One such phenomenon is thermal bowing of 
the metallic TPS panels. 3 Bowing can occur because 
metals have relatively large coefficients of 
thermal expansion and because in the course of 
performing their function as a thermal protection 
system the TPS panels must sustain large tempera- 
ture gradients through their thickness. These 
conditions lead to large thermal expansion of the 
outer portion of the panel and small thermal expan- 
sion of the inner portion of the panel. The square 
panels are held down at the corners but otherwise 
thermal expansion is unrestrained. As a result the 
center of the panel bows out into the stream so 
that the vehicle moldline is altered. An array of 
panels would take on a quilted configuration. 

Analytical and experimental studies of two- 
dimensional flow over a wavy»wall are available in 
the literature, however, extrapolation of these 
results to the three-dimensional flow over an array 
of bowed panels or spherical dome protuberances 
results In an Inaccurate description of the 
phenomena. The analysis must be done in three- 
dimensions to adequately predict the physical 
flows. No analytical studies of a configuration 
resembling the thermally bowed TPS were found in 
the literature, probably because large amounts of 
computer time and storage are required to perform 
the analysis. Consequently, this preliminary study 
was undertaken to determine if the Increased 
pressure and heating effects were significant and 
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if the problem required in-depth consideration. 
Thermally bowed panels were modeled by the generic 
configuration of spherical domes. Numerical grid- 
ding and computer storage limitations prohibited 
modeling an array of full scale TPS panels so a 
single row of 4.6 inch diameter domes was selected 
to be conpatlble with the grlddlng requirements. 
Dome heights of 0.1, 0.2, and 0.4' Inches were 
considered In a flow field with a laminar-boundary- 
layer thickness of 0.2 Inches. Free stream 
conditions are for Mach 7 flight at 120,000 ft. 
altitude. The analysis of this flow configuration 
will yield dimensionally similar solutions. 

However, the proper balancing of the similarity 
parameters to ascertain that the small scale pro- 
cesses in the boundary layer will fall Into the 
similarity pattern requires additional analytical 
and experimental data not yet available . 4 
Therefore, no extension of the results of this 
study to larger protuberances will be made at this 
time. 


The analysis employs a three-dimensional 
formulation of the time dependent equations for 
conservation of mass, momentum, and energy. A 
MacCormack explicit time-split predictor-corrector 
finite difference algorithm was used to obtain a 
steady-state solution on Langley's Cyber 203 vector 
processing computer. The results of those solu- 
tions are presented In this paper using state-of- 
the-art computer graphics techniques. 

A comprehensive experimental program to verify 
the results obtained In this analysis Is scheduled 
for a later time. Some preliminary small scale 
tests have provided qualitative correlation with 
some of the predicted phenomena; no experimental 
data will be presented in this paper. 


Analysis 


Governing Equations 

Aerothermal loads on spherical domed protuber- 
ances were computed by' mathematically modeling the 
flow field with the continuum mechanics equations 
of motion. Air, the fluid medium, was modeled as a 
compressible, viscous, thermally conducting gas. 
Limiting assumptions Include considering air as an 
isotropic, newtonlan ideal gas in local thermo- 
dynamic equilibrium. The three-dimensional 
time-dependant system of equations describing this 
model include the continuity equation (conservation 
of mass), the Navler-Stokes equations of motion 
(conservation of momentum), the energy equation 
(conservation of energy), and the Ideal gas 
equation (thermodynamic state equation). These 
equations have been derived in the literature 3 and 
are shown below written In cartesian tensor form 
where repeated Indices Indicate a summation over 

3(pwJ 

the range of Indices (l.e. , — gjj — - for 1 » 1, 2, 3 
3pW, 8pW ? 8pW r 1 

Is In addition w lt w 2 , W 3 

12 3 

are the velocity components u, v, w, and Xj, X 2 , X 3 
are the coordinate directions x, y, z. The 
Kroneckeh delta, equals 1 when 1 - j and 

equals zero when 1 f j. 
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(results ifji three equations, one for each 
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Energy equation . 
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for 1, j, k, ■ 1, 2, 3 


Ideal gas state equation . 

P = pR g T (4) 

These equations along with the ideal gas 
relationship, 

e = c v T (5) 

and the Sutherland viscosity formula, 

T 3/2 .. ,„-8 lb sec 

" 2,270 T + I5S7F x 10 f 2 

constitute a set of five partial differential 
equations and one algebraic equation in six 
dependant variables: three velocity components, u, 

V, and w; density, p; temperature, T; and pressure, 
P. 


The physical domain of the problem is shown in 
figure 1* The thermally bowed TPS is modeled as a 
series of doroes on a flat plate. Because of 
symmetry planes at the center of each dome and 
between domes, only one-half of a dome need be 
modeled. The dome was specified by its height, h, 
above the flat plate and its radius, r 0> at the 
dome/flat plate Intersection. The Intersection 
between the dome and flat plate was made 
mathematically smooth by modeling a fillet region 
at r 0 ± e with a Hermite polynomial. 

Boundary Conditions 

A boundary layer profile was input as the 
upstream boundary plane data and held constant 
throughout the computations. The boundary layer 
profile was determined using a two-dimensional 


boundary layer analysis 6 for flow over a flat 
plate. The flow configuration is shown in figure 1 
and free stream conditions for Mach 7 flight at an 
altitude of 120,000 ft. are listed in Table I, The 
boundary layer profiles for the starting plane are 
shown in figure 2. These profiles, along with a 
crossflow velocity, w, of zero were impressed on 
the starting plane (x = 0) and held constant 
throughout the computations. For a flat plate with 
no dome the boundary layer height, fif p » at the^ 
center of the three-dimensional region (x “ x c J 
0.2 inches. Dome heights, h, with ratios h/o fp = 

0. b, I, and 2 were investigated. The ratio of dome 
radius to boundary layer height, r 0 /6f p , was held 
constant at 11.5. 

Doth side boundaries are planes of symmetry on 
which the dependent variables are determined by 
quadratic extrapolation from interior points so 
that the cross flow gradients, 0/9z, of all state 
variables are zero and the cross flow velocity, w, 
is zero. The upper boundary and downstream 
boundaries are specified as no change boundaries, 
i.o., set equal to the plane next to them. The 
flat -pi ate/ dome surface is a no-slip boundary with 
a constant wall temperature of 1440°R. This wall 
temperature is both the upper use temperature and 
the flat plate radiation equilibrium temperature 
for the input boundary conditions for a titanium 
TPS assuming a high emittance coating and zero 
thermal diffusivity. 

Computational Domain 

In order to generate a boundary fitted 
coordinate system and provide increased resolution 
in the boundary layer, the physical domain (fig. 1) 
was transformed into a rectangular parallelepiped 
computational domain using a two-boundary grid 
generation technique. The computational domain, 
shown in figure 3, has equally spaced grid lines in 
the C* n, and r, directions. 

Transformation of the governing differential 
equations, eqs. 1, 2, and 3, to the computational 
domain requires replacement of all spatial 
derivatives, 9/OX.j for i = 1, 2, 3, with the 
following. 


* i nc „ 

a a i 

nxT = 5cT SxT 


(7) 


where C £ for £ = 1, 2, 
coordinates C, n and C, 

0 3 oc a Qn_ 

ax. “ oc ax. On ox. 

i i i 


3 are the transformed 
i .e. , 

0? The resultin 9 
‘ i 


equations of motion are* 
Continui ty, 


HE 

Ot 


+ 


0 (pw. ) oC 
~ac W~ 


£ 1 



Navier Stokes equations of motion 
a(pw.) 0 (pw.w.) OC 9 (t..) OC 


( 9 ) 


ORIGINAL PAGE US 
OF POOR QUALITY 


where 
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Written in short form the system of partial differ- 
ential equations Is, 
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The required Jacobian transform matrix. 
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is determined by a linear mapping of the boundaries 
at n= 0 and n= 1 of the computation domain onto 
the boundaries of the physical domain at y = 0 and 
y = ymax In the following form; 


Y = Y (5, 1, S) n + Y (5. 0, c) (1 - n). (13) 


3 


ORlGfttftL PA.GS S» 

OF POOR QUALITY 


further transformation to Increase resolution In 
the boundary layer region entails an exponential 
grid stretching In the n direction of the form. 


e kn -l 
e k -1 


(14) 


where k is a free parameter that adjusts grid 
spacing. Equation 13 becomes, 

V « Y (€, 1, C) n + Y (5, 0, C) (1 - n) (15) 

For the geometry of the subject problem the 
resulting Jacobian matrix Is of the form, 
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where §, -g*. and •§, which vary for the flat 

plate, fillet, and dome regions as defined by the 
radius measured from the dome centerline. 
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are shown in the Appendix. 

These Jacobian elements are calculated and 
stored for each grid point in the domain thus 
allowing the solutions from the computational 
domain to be transformed to the physical domain. 

Substituting the transform Jacobian (ea. 16) 
into the short form of the equation system (eq. 11) 
reduces the system to the following, 

.aU . at 1 + 5E.la + a&5n + i!i3n + ^i —L. = o 
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Numerical Procedure 

A finite-difference solution to the governing 
equations was attained using the second-order 
accurate MacCormack explicit time-split 
predictor-corrector algorithm. The algorithm 
splits the differencing scheme into a series of 
one-dimensional operators. The operator in the 
principle flow direction, Lp, advances the 
computation by one time increment while the 
operators in the other two coordinate directions, 

L and L c , are divided into two steps of one-half a 
tPme-increment each and applied symmetrically 
around the principle flow direction operator. This 
arrangement allows the principle flow direction 
operator to run at a Courant, Friedricks, Lewy 
(CFL) number close to the optimum value of one. 

The operators are applied serially as follows: 


Each operator has two steps, a predictor step that 
advances the solution by its time increment based 
on a backward spatial differencing of the 
one-dimensional equations and a corrector step that 
recalculates the advancement based on a forward 
spatial differencing of the predicted results and 
averages the result with the predictor. The 
predictor must therefore lag the corrector by a 
spatial increment in each coordinate direction. 

The combination of these two steps results in a 
second-order-accurate differencing scheme in time 
and space. As an example the first operation 
LjUt/2) applied to the data at time step Ms as 
follows, 


Predictor step: 
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Corrector step : 
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Where bars over the variables indicate predictor 
values and asterisks over the variables indicate an 
intermediate result in the serial application of 
operators as shown in equation 19. The other 
operators have similar form. 

Although a complete stability analysis of this 
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where c is the local speed of sound. 

In the formulation of this algorithm terms of 
order three and higher have been truncated. These 
terms are not significant in many flow regimes, 
but, when modeling high energy flow fields with 
strong shocks, as in this paper, they can be 
significant. In the region of a shock, where 
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